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A hybrid-exchange density-functional theory study of the electronic structure of 
MnV 204 : Exotic orbital ordering in the cubic structure 
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The electronic structures of the cubic and tetragonal MnV 2 04 have been studied by using hybrid- 
exchange density functional theory. The computed electronic structure of the tetragonal phase shows 
an anti-ferro orbital ordering on V sites and a ferrimagnetic ground state (the spins on V and Mn are 
anti-aligned). These results are in a good agreement with the previous theoretical result obtained 
from the local-density approximation+f/ methods [S. Sarkar, et. ah, Phys. Rev. Lett. 102, 216405 
(2009)]. Moreover, the electronic structure, especially the projected density of states of the cubic 
phase has been predicted with a good agreement with the recent soft x-ray spectroscopy experiment. 

Similar to the tetragonal phase, the spins on V and Mn in the cubic structure favour a ferrimagnetic 
configuration. Most interesting is that the computed charge densities of the spin-carrying orbitals 
on V in the cubic phase show an exotic orbital ordering, i.e., a ferro-orbital ordering along [110] but 
an anti-ferro-orbital ordering along [TlO]. 

PACS numbers: 71.15.Mb, 75.47.Lx, 75.50.Gg, 75.25.Dk 


I. INTRODUCTION 

Many fascinating phenomena in condensed-matter 
physics were discovered in transition-metal oxides 
(TMOs) and their closely related compounds such as 
the first high transition temperature superconductivity 
in YBaCuOs Eli- An important ordering phenomenon, 
Orbital ordering (00) has a long history in the physics 
of TMOs, back to the 1930s, when the concept of charge 
ordering in Fe 304 was proposed Q. 00 occurs when the 
orbital degeneracy is lifted essentially by the interaction 
with the lattice environment Bi- The best-known ex¬ 
ample of 00 is the Jahn-Teller (JT) distortion observed 
in LaMnOs Q. Therein the eg degenerate manifold 
and dx 2 _y 2 ) in a cubic environment is further splitted by 
the JT distortion owing to partially filling of electrons. 
The underlying physics in 00 has been well accounted 
for by the so-called Kugel-Khomskii model Q • Most ex¬ 
perimental and theoretical works so far have been fo¬ 
cused on the system with partially filled eg orbitals, i.e., 
eg active. This type of orbital occupancy can lead to 
the strongest JT effect because the related d-orbitals di¬ 
rectly point to ligands and strongly hybridise with the 
2p-orbitals of the anion. 

Recently much attention has been paid to vanadium 
spinels AB2O4 (A = Mg, Co, Fe, Mn, and Cd, etc and 
B = V) |§-[l^ owing to their fascinating spin and orbital 
orderings accompanied by complex structural transitions 
at low temperature. The interplay between structure, or¬ 
bital, and spin implies that the underlying physics could 
be very interesting, and that there might be large applica¬ 
tion potential in these materials. For example, artificially 
tuning structure by applying external magnetic field can 
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be thought of as a prototype of TMO-based metamate¬ 
rials or even quantum metamaterials [HI [3. Another 
unique point for vanadium spinels is that they are t 2 g- 
active, in sharp contrast to most of the known TMOs 
showing 00. In this type of compound, the ions on 
the B sites form a pyrochlore lattice. The two d-electrons 
on V occupying two t 2 g orbitals lead to a total spin of 
S = 1, resulting in a t 2 ^-active 00. The 00 in this type 
of compounds is further complicated by their intrinsic 
geometrical frustration, which played an important role 
in 00. In 1939, the charge ordering on the B-site py¬ 
rochlore lattice was proposed to cause a sharp increase 
of the electrical resistivity when cooling below ^ 120 K 
in Fe 304 (the ’Verwey’ transition) Q. Anderson was 
the first to realise the intimate relationship between the 
proton ordering in the ice, the Verwey charge ordering, 
and the ordering of Ising spins M- Regardless of the 
specific material types, geometrical frustration plays a 
crucial role to determine the spin or charge configura¬ 
tions, by minimising exchange or Coulomb energies. This 
will in turn trigger many interestin g ph enomena such as 
orbital-glass and orbital-ice states |16l-[lq. In addition, 
the much faster response of electron charges as compared 
to electron spins will find many potential _^plications in 
advanced electronics such as orbitronics |1^ . 

MnV 2 04 (See FigH]), in which the Mn^^ ion is in a d^ 
half-filled high-spin configuration {S = |), first experi¬ 
ences a magnetic transition to the collinear ferrimagnetic 
state at T = 56 K, and then a structural distortion at a 
slightly lower temperature T = 53 K, along with a tran¬ 
sition to the non-collinear ferrimagnetic state [lO, 1^ . 
The observed magnetic ordering below T = 56 K is ferri¬ 
magnetic, i.e., the magnetic moments on Mn^^ and 
are anti-aligned. In the structural transition, the sym¬ 
metry is lowered from cubic to tetragonal. The 00 of 
the ground state in the tetragonal MnV 204 structure 
has been shown theoretically as an A-type antiferro- 




2 


orbital ordering with a propagation vector of (002). 
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FIG. 1: (Color online.) the conventional cells in the cubic 
(left) and tetragonal (right) MnV 2 04 crystal structures are 
shown. Mn is depicted as small yellow ball, V as large green 
ball, and O as large red ball. 


Previously, the electronic structure and magnetic prop¬ 
erties of the tetragonal MnV 2 04 have been investigated 
theoretically using LDA (local-density approximation) + 
U method [l^, where an orbital ordering among was 
proposed, along with a non-collinear magnetic ordering. 
Following this first-principles calculation, an analytical 
modelling [22| has further confirmed that an antiferro- 
orbital ordering exists in the tetragonal MnV 204 . On 
the other hand, the electronic structure (especially 00) 
and magnetic properties of the cubic MnV 204 are yet 
to be studied in detail. It might be worthwhile to (i) 
employ a computational method without any adjustable 
parameters, (ii) take into account the electron correlation 
properly and (hi) compare the electronic structures of the 
cubic and tetragonal MnV 204 . Hybrid-exchange func¬ 
tional theory (HDFT), in which the exact exchange is hy¬ 
bridised with generalised gradient approximation (GGA) 
functional to localise d-electrons, has performed well for a 
wide range of inorganic and organic compounds mil. 
In this paper, HDFT with PBEO functional has been 
used to study the electronic structure and magnetic prop¬ 
erties of the cubic and tetragonal MnV 204 . The results 
presented here are not only in a good agreement with the 
previous theoretical and recent experimental studies, but 
also point to an exotic 00 state that mixes ferro-orbital¬ 
and anti-ferro-orbital orderings. The rest of the discus¬ 
sion is organised as the following: in ^IIl computational 
details are introduced, in m the calculation results are 
presented and discussed, and in m some general con¬ 
clusions are drawn. 


II. COMPUTATIONAL DETAILS 

The calculations of the electronic structures of the 
tetragonal and cubic MnV 204 were carried out by us¬ 


ing DFT and hybrid-exchange functional PBEO as im¬ 
plemented in the CRYSTAL 09 code [2^. The crystal 
structures of the cubic (space group Fd3m) and tetrago¬ 
nal (space group Mi/amd) MnV 204 experimentally de¬ 
termined in Ref. ini have been adopted here to perform 
all the calculations. The basis sets for the atomic or¬ 
bitals centred on the Mnlai, VISH , and 0|33 atoms, 
which are designed for solid-state compounds, were used. 
The Monkhorst-Pack samplings [33| of reciprocal space 
were carried out choosing a grid of shrinking factor to 
be 6 X 6 X 6 (6x6x5) in order for a consistency 
with the ratios among reciprocal lattice parameters for 
the cubic (tetragonal) MnV 2 04 . The truncation of the 
Coulomb and exchange series in direct space was con¬ 
trolled by setting the Gaussian overlap tolerance crite- 
ria to 10-®,10-®,10-®,10-®, and The self- 

consistent field (SCF) procedure was converged to a tol¬ 
erance of 10“^ a.u. per unit cell (p.u.c). To accelerate 
convergence of the SCF process, all the calculations have 
been converged by using a linear mixing of Fock matrices 
by 30%. 

Electronic exchange and correlation are described us¬ 
ing the PBEO hybrid functional [25|, free of any empir¬ 
ical or adjustable parameter. The advantages of PBEO 
include a partial elimination of the self-interaction er¬ 
ror and balancing the tendencies to delocalize and lo¬ 
calize wave-functions by mixing a quarter of Eock ex¬ 
change with that from a generalized ^adient approxi¬ 
mation (GGA) exchange functional [^. The broken- 
symmetry method [sj is used to localize collinear op¬ 
posite electron spins on atoms in order to describe the 
anti-ferromagnetic state. 


III. RESULTS AND DISCUSSIONS 
A. Projected densities of states 

The projected densities of states (PDOS) of the cubic 
and tetragonal MnV 2 04 structures for the AEM config¬ 
uration (the spins on Mn and V are anti-ligned) have 
been shown in Eigl2l The zero energy is aligned with the 
valence band maximum (VBM). 

In the cubic and tetragonal structures, Mn 3d PDOS 
(Eigl2^ and b) shows that the five d-orbitals including 
dx 2 _y 2 ^ dz 2 , dxz^ dyz^ and dxy are singly occupied, thus 
giving a spin-1. Erom the analysis of Mulliken spin den¬ 
sities and V 3d PDOS (EiglSJ: and d), the two electrons 
occupying t 2 g states can be identified as the origin of the 
00 both in the cubic and tetragonal MnV 204 . In the cu¬ 
bic phase, the molecular orbitals delocalized among the 
three t 2 g states are occupied, whereas in the tetragonal 
phase, only dxz and dyz are involved. The O 2p PDOS 
(EigISt) for the cubic phase are much more delocalized 
than those for the tetragonal (Eigl2f); this might be due 
to the elongation of the lattice vector a and b in the 
tetragonal structure. The O 2p PDOS are particularly 
dominant at ^ 5 eV below the VBM for both the cubic 





















3 


Cubic 


Tetragonal 


Cubic MnV^O^ AFM Mn-3d-orbital PDOS Tetragonal MnV 204 AFM Mn-3d-orbital PDOS 



Energy (eV) Energy (eV) 


(a) 


(b) 


Cubic AFM V-3d-orbital PDOS 


Tetragonal MnV 204 AFM V-3d-orbital PDOS 



-10 -5 0 5 10 

Energy (eV) 



-10 -5 0 5 10 


Energy (eV) 


(c) 


(d) 


Cubic MnV 204 AFM 0-2p-orbital PDOS 


Tetragonal MnV 204 0-2p AFM PDOS 



Cubic MnV 204 AFM PDOS 


Tetragonal MnV 204 AFM PDOS 



Energy (eV) Energy (eV) 

(g) (h) 


FIG. 2: (Color online.) The PDOS of the 3(i-orbitals on Mn, the 3(i-orbitals on V, and the 2]2-orbitals on O of the cubic (left) 
and tetragonal (right) MnV 204 are shown. The PDOS onto d^ 2 _y 2 is in green, in red, dxz in blue, dyz in purple, and dxy 
in black, respectively. For 2p-orbitals, the PDOS onto 2px is in green, 2py in red, and 2pz in black, respectively. In (g) and 
(h), the total PDOS onto Mn d-orbitals is in red, V d-orbitals in green, and O 2p-orbitals in black, respectively. 
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and tetragonal structures, which overlap with the rather 
weak Mn 3d and V Sd PDOS. The O 2p orbitals involved 
here will give rise to the so-called oxygen bonding states 
that feature the hybridisation between the d-orbitals on 
transition metals and p-orbitals on O atoms. In addition, 
it can be observed by comparing the spin-up and spin- 
down PDOS, that the O 2p states have been strongly 
spin-polarized by the magnetic moments on Mn and V. 
The PDOS onto the V sites have also been compared 
with the previous results obtained by the LDA+I/ HI, 
which suggests that there is a qualitative agreement be¬ 
tween them except that the DPT band-gap computed 
here 3 eV) is larger than the previously computed 
one 1.1 eV) 0 and the observed 1.1 eV)®. 
The computed O 2p PDOS are in good agreement with 
the intensities measured in the recent K-edge x-ray ab¬ 
sorption and emission spectra; the combination of the 
theoretical work and experiments will be published in a 
forthcoming paper j^ . 

The computed crystal-field splitting between t 2 g and eg 
states of Mn^^ and in the cubic MnV 204 can be read 
from Fig 12^, which is ^ 1 eV. On the other hand, in the 
tetragonal MnV 2 04 , for both Mn^^ and dj. 2 _y 2 and 
dxy are closely aligned while the other three Sd-orbitals 
overlap well, as shown in Figl2]3. This picture is consis¬ 
tent with the previous calculations reported in Ref. 0 . 
The crystal-field splitting between these two groups is 
^ 1 eV. The on-site Coulomb interaction U can be ap¬ 
proximated by observing the gap between lower Hubbard 
d-bands and uppers ones 5 eV for Mn site, and ^ 7 
eV for V), which is much larger than the crystal-field 
splitting computed here. This will result in a high-spin 
state for Mn^^ and ions, which is consistent with 
the previous experiments 0 . 


B. Exchange interactions 

The magnetic structure could be much more compli¬ 
cated owing to the exchange interaction between the 
nearest-neighbouring Mn (V) atoms 0 and the spin- 
orbit coupling (SOC), which will cause the non-collinear 
magnetic ordering. However, this topic is beyond the 
scope of this paper that is focused on the electronic 
structure and the collinear magnetic ground state. The 
Mulliken spin densities on the Mn and V sites for 
the cubic (tetragonal) MnV 204 are 4.8 (4.7)/iB and 
—2.0 (— 1.9)/iB5 respectively. They are close to the ex¬ 
pected values, i.e. bjiB for Mn and —2pB for V (anti¬ 
aligned). For the cubic (tetragonal) structure, the com¬ 
puted total energies for a conventional cell with the spins 
on Mn and V aligned (ferromagnetic) is higher than the 
anti-aligned (AFM) by 937 (754) meV. This in turn gives 
an approximate exchange interaction of ^ 38 meV (di¬ 
vided by the factor of 4 x S'Mn x 5'v=20, where 4 is the 
number of Mn-V pairs in the unit cell for the tetragonal 
structure). Similarly we can estimate the exchange in¬ 
teraction for the cubic structure, which is ^ 23 meV. On 


the other hand, the exchange interaction between the Mn 
and V spins can also be estimated by using the super¬ 
exchange formalism, where t can be quantified in the 
PDOS, which is ^ 0.5 eV for Mn and ^ 1 eV for V. By 
using these hopping integrals, the exchange interaction 
can be estimated as {JJ = which is ^ 59 

meV. This is in the same order as that calculated from 
the total energy differences aforementioned (here an av¬ 
erage value is taken for on-site Coulomb interaction). A 
detailed discussion about the exchange interaction, tak¬ 
ing into account all the d-orbitals on both Mn and V 
would be great, but beyond the scope of this paper. 


C. Orbital ordering 

The 00 in the cubic and tetragonal MnV 204 can be 
illustrated by the spin densities on V (see Figl3])- the dif¬ 
ference between the charge densities of spin-up and spin- 
down orbitals (essentially the charge densities of the spin- 
carry orbitals). In the tetragonal structure, the relative 
orientation rotation of neighbouring orbitals is illustrated 
by the red arrows in FiglSb; this is consistent with the 
previous calculations presented in Ref.[T2| in which an 
anti-ferro-orbital ordering (AFOO) has been predicted. 
The orbital orientations for the nearest-neighbouring 
(NN) orbitals are perpendicular to each other (defined 
as AFOO). In sharp contrast, for the cubic structure the 
NN orbtibals (labeled by X and Y in FiglS^t) are orga¬ 
nized in an exotic way, in which the orbitals on one of 
the four symmetry-inequivalent V atoms in the unit cell 
has a different orbital orientation (labeled by Y), perpen¬ 
dicular to the others (labeled by X). The orbital orien¬ 
tations corresponding to Figl3^ are further illustrated in 
FigISji, e, and f, which are the views of the spin densities 
from the lattice direction [100] (a), [010] (6), and [001] 
(^, respectively. This leads to an exotic 00: a ferro- 
00 (FOO) along [110] (blue arrow in Figl3]), whereas an 
AFOO along [110] (red arrow in FigHJ. This 00 can 
probably be attributed to (i) the crystal field environ¬ 
ment formed by neighbouring oxygen atoms and (ii) the 
Coulomb interaction between d-orbitals. This interesting 
00 has been further confirmed by the calculations with a 
different type of exchange-correlation functional, B3LYP 
j35j , which relies on the empirical parameters for mixing 
exact exchange and GGA exchange functional. The spin 
densities on V predicted by using B3LYP functional are 
shown in Figl3t, in which FOO is along [110] and AFOO 
[110] (all the atoms are made partially transparent in 
order to illustrate 00). 

To see the orbitals involved in the 00 of the cubic 
MnV 2 04 more clearly, the atomic orbital compositions of 
the bands (see FigH]) at the F-point that contribute most 
to the two peaks at —0.26 eV and —0.64 eV in the V spin- 
down PDOS (see FiglJt) have been investigated. The 
bands VB—8 and VB—9 contribute the most to the peak 
at —0.26 eV, whereas VB—19 and VB—20 to the one at 
—0.64 eV. Although the atomic coefficients for d-orbitals 
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FIG. 3: (Colour on line.) The spin densities on V of the AFM state in the conventional unit cell of the cubic (a) and 
tetragonal (b) structures of MnV 2 04 are shown. In the tetragonal an anti-ferro-orbital ordering is found. In contrast, in the 
cubic structure, the rotation of the orbital orientation illustrates an exotic 00, in which an ferro-orbital ordering is along [110], 
while anti-ferro-orbital ordering along [110]. This exotic 00 is further confirmed by the calculation from B3LYP functional; 
the resulting spin densities are shown in (c). The B3LYP spin-density calculation is performed in a unit cell, but shown in a 
conventional cell. The views from [100] (a), [010] (5), [001] lattice directions are shown in (d), (e), and (f), respectively. 
The isovalue is chosen as 0.04e/A^. 


are slightly different among V atoms, approximately the 
most dominant linear combinations of d-orbitals for each 
V atoms contained in these bands read, 



FIG. 4: (Golour on fine.) The spin-down band structure 
(from -0.8 to 0 eV) of the AFM configuration in the cubic 
MnV204 is shown. Note that the two red arrows refer to 
VB-8 and VB-9, and VB-19 and VB-20, respectively. 


0VB-8 = ^•‘^\dxz) + ^-Mdyz) + 0-2|<iccy), (1) 

0VB-9 = -0.l\dxz) + 0.2\dyz) - 0.2\dxy), (2) 

0VB-19 = 0.2\dxz) - 0.1\dyz) - 0.2\dxy), (3) 

0VB-2O = ^‘‘^\dxz) + 0.2|(i^^) + 0.1\dxy). (4) 

The corresponding spherical harmonic functions have 
been plotted in Figl5] to see the contributions of d- 
orbitals. From the perspectives of linear combination 
of atomic orbitals (LCAO), the strong variations of the 
orbital orientation illustrated in Figl3] and FiglSl result 
from the minimisation of the total energy respective to 
atomic orbital coefficients. On the other hand, the pro¬ 
cess of energy-minisation can be seen as the orbital re¬ 
orientation driven by Coulomb and crystal-field interac¬ 
tions. 

There are two concepts that are closely related to this 
exotic orbital ordering, including orbital glass and ice 
states. The orbital glass state is formed by the randomly 
orientated orbitals, which result from Jahn-Teller distor¬ 
tion and orbital frustrations in a similar manner to spin 
frustrations [iSlil- In light that orbital glass states 
have been observed in the spinel compounds FeCr 2 S 4 
and LuFe 204 [mil (which has similar chemical for¬ 
mula to a spinel, but a different structure), this is not 
completely unexpectable in the cubic MnV 204 . There¬ 
fore, these calculations could suggest that at T = 53 
K, there is not only a structural transition, but also an 
orbital-ordering transition from the exotic 00 illustrated 
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FIG. 5: (Colour online.) The linear combination of harmonic spherical functions are shown. 2Yxz -\-Yyz + ‘2Yxy (corresponding 
to eq.([T])) are shown in (a), —Yxz + ^Yyz — 2Yxy (corresponding to eq.(|2])) in (b), 2Yxz — Yyz — 2Yxy (corresponding to eq.(|3])) in 
(c), and 2Yxz + ‘2.Yyz + Yxy (corresponding to eq.JH)) in (d). Notice that the linear combination illustrated in (c) changes the 
orientation of orbitals significantly away from their original ones, which is closely related to the exotic 00 shown in Fig|3^. 


in Figl3]to anti-ferro orbital ordering [a , when lowering 
temperature. In addition, Chern and the collaborators 
have presented an ice model consisting of triplet orbital 
variables (such as p-orbitals), which is however closely 
related to orbital-driven many-body phenomena in opti¬ 
cal lattice [TsI. The analogue of spin-interacting Hamil¬ 
tonian, an orbital-exchange model describing Coulomb 
interactions has been used therein. The exotic ordering 
presented here is closely related to the newly proposed 
concept of the orbital ice 0 , similar to the spin ice. 

IV. CONCLUSION 

The electronic structures of the cubic and tetrago¬ 
nal MnV 204 structures have been computed by using 
hybrid-exchange density functional theory PBEO and 
B3LYP. The calculated PDOS of the tetragonal struc¬ 
ture has been in a good agreement with previous theoret¬ 
ical results. The most important finding here is that the 
charge densities of spin-carrying orbitals suggest a pos¬ 
sible existence of an exotic orbital ordering: FOO along 
[110] and AFOO [110] in the cubic MnV 204 . The theo¬ 
retical prediction presented here could be validated in the 
future experiment. For example, the different 00s along 
those two lattice orientations aforementioned could lead 


to different electron transport properties. More advanced 
theoretical methods such as dynamical mean-field theory 
[38| is also needed to further understand the electronic 
structure of MnV 204 . In the ground state the spins on 
Mn and V are predicted to be anti-aligned, suggesting a 
ferrimagnetic state for both structures, which is in good 
agreement with previous theoretical and experimental re¬ 
sults. The lower and upper Hubbard bands of Mn and 
V d-electrons have been shown clearly below and above 
the band gap. The on-site Coulomb interaction can be 
estimated by using the gap between the lower and upper 
Hubbard bands, which is ^ 7 eV for V and ^ 5 eV for 
Mn, respectively. The transfer integral responsible for 
the delocalization is in the order of 1 eV, which should 
be attributed to a combination of the first-order direct 
hopping between Mn and V and the second-order effect 
via O atoms. 
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